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Type la supernovae, the thermonuclear explosions of white dwarf stars composed of 
carbon and oxygen, were instrumental as distance indicators in establishing the accel- 
eration of the universe's expansion. However, the physics of the explosion are debated. 
Here we report a systematic spectral analysis of a large sample of well observed type la 
supernovae. Mapping the velocity distribution of the main products of nuclear burn- 
ing, we constrain theoretical scenarios. We find that all supernovae have low-velocity 
cores of stable iron-group elements. Outside this core, nickel-56 dominates the super- 
nova ejecta. The outer extent of the iron-group material depends on the amount of 
nickel-56 and coincides with the inner extent of silicon, the principal product of in- 
complete burning. The outer extent of the bulk of silicon is similar in all SNe, having 
an expansion velocity of ~ 11000 km s"^ and corresponding to a mass of slightly over 
one solar mass. This indicates that all the supernovae considered here burned similar 
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masses, and suggests that their progenitors had the same mass. Synthetic light curve 
parameters and three-dimensional explosion simulations support this interpretation. 
A single explosion scenario, possibly a delayed detonation, may thus explain most type 
la supernovae. 

When a white dwarf (WD) composed of carbon and oxygen accreting mass from a companion star in 
a binary system approaches the Chandrasekhar mass [Mcb ~ 1-38 solar masses (M©)], high temperature 
causes the ignition of explosive nuclear burning reactions that process stellar material and produce en- 
ergy. The star explodes leaving no remnant, producing a Type la supernova (SNIa) (i ). At high stellar 
material densities, burning yields nuclear statistical equilibrium (NSE) isotopes, in particular radioactive 
^^Ni which decays to ^^Co and ^^Fe making the SN bright (2). At lower densities intermediate mass el- 
ements (IME) are synthesised. Both groups of elements are observed in the optical spectra of SNela Q. 
An empirical relation between an observed quantity, the B-magnitude decline over the first 15 days after 
maximum luminosity [Anii5(B)] and a physical quantity, the SN maximum luminosity (4"), can be used to 
determine the distance to a SNIa. This method was applied to very distant SNela leading to the discovery 
of the accelerating Universe (Old]). How the explosion actually proceeds is however debated (^|l]|9ll, as is 
the nature of the progenitor system: Accretion may occur either from a more massive companion (e.g. a 
giant) or via the merging of two carbon-oxygen WDs (70). This casts a shadow on the reliability of SNela 
as distance indicators, as intrinsically very different explosions may result in the observed correlation. 

We derive the distribution of the principal elements in 23 nearby SNela (distances <40Mpc) with 
good spectral coverage extending from before maximum to the late nebular phase, about one year later. 
The sample (table SI) covers a wide range of light-curve decline rates and includes peculiar objects like 
SN2000cx, which violates the luminosity-decline-rate relation. 

Because of the hydrodynamic properties of the explosion, the expansion velocity of the ejecta is pro- 
portional to radius and serves as a radial coordinate. As the SN expands, deeper layers are exposed. The 
outer layers, visible in the first few weeks after the explosion, are dominated by IME. Because silicon 
is the most abundant IME, we measured its characteristic expansion velocity from the blueshift of the 
absorption core of the strong Sill 6355 A line in all spectra where it was visible. This velocity decreases 
with time. Fitting the postmaximum velocity evolution and extrapolating it to the earliest times, when the 
outermost parts of the ejecta are visible, we derived the outer extension of the bulk of Si. This represents 
a lower limit of the outer extent of burning. The Si velocity !;(Si) is similar in all SNela, regardless of 
their luminosity: v{Si) = 11900 ± 1300km s~'. 

A few SNe, defined as High Velocity Gradient (HVG) SNe ((77|) . are responsible for most of the disper- 
sion. They have a rapidly decreasing i;(Si) before maximum and very high-velocity Call lines, possibly 
the result of high- velocity blobs that carry little mass and kinetic energy but cause High- Velocity absorp- 
tion Features (HVF, (I72]| ) that can abnormally broaden the Sill line profile il2\\13h Although excluding 
them from the sample decreases the dispersion in v{Si) significantly [t;(Si) = 11300 ± 650km s"^], we 
include all these SNe in our discussion. 

The maximum Si velocity thus measured is a conservative estimate. The deep absorption core is 
produced in layers of high silicon abundance. Silicon is present at higher velocities, indicated by the 
wavelength of bluest absorption in the earliest spectrum of each SN. However, measuring the bluest 
absorption velocity yields a large scatter because the earliest spectra have different epochs and may be 
affected by HVFs {12 ). Our method reliably determines the outer location of the bulk of IMEs. 

The inner extent of Si, determined from the asymptotic velocity of the Sill 6355 A Une in post- 
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maximum spectra [fig. 1 in (I77])1. is a steep function of Ami5(B). The brightest (slowest declining) 
SNe have the thinnest Si zones. 

The inner ejecta, dominated by NSE elements, are best observed ~1 year after the explosion, when 
dilution caused by expansion makes the SN behave like a nebula, exposing the deepest layers. Collisions 
with the fast particles produced by the decay ^^Ni — > ^^Co ^^Fe heat the gas, which cools emitting 
radiation mostly in forbidden lines. 

We modeled the nebular spectra using a code that computes line emission balancing heating and 
cooling in non-local thermodynamic equilibrium (NLTE) (14), including density and abundance strat- 
ification. We adopted the density-velocity distribution of the standard, one-dimensional, Mch explosion 
model W7 (7), (fig. SI). In the nebular phase the gas is transparent, and line emissivity depends on the 
mass of the emitting ion. Accurately estimating this mass requires determining the ionization state of the 
gas. Forbidden lines of Fell and Felll dominate SNIa nebular spectra, reflecting the high abundance of 
NSE material. Fe is mostly the product of ^^Ni decay, which provides heating. The stable, neutron-rich 
isotopes ^^Fe and ^^Ni do not contribute to heating, but do contribute to cooling, because they also emit 
forbidden lines. Their presence affects the ionization balance. Both are mostly produced deep in the 
WD, at the highest densities. The ^'^Fe nebular lines have wavelengths indistinguishable from those of 
5^Fe. 

We determined the mass and distribution in velocity of the Fe isotopes (and thus of ^^Ni) simultane- 
ously fitting the ratio of the two strongest Fe emissions(fig. SI). One, near 4700 A includes both Fell 
and Felll lines, while the other, near 5200 A, is only due to Fell. A low upper limit to the mass of ^^Ni is 
set by the absence of strong emission lines (in particular at 7380 A). 

In the deepest layers (Fig. [T]), all SNe contain ~0.1 to 0.3 M© of stable NSE isotopes, with a large 
scatter and no dependence on Ami5(B) [see also U5i)]. 

As expected, the ^^Ni mass correlates inversely with Ami5(B), ranging from 0.9 M© for the slowest 
declining (most luminous) SNe to 0.1 Mq for the fastest declining (dimmest) ones. The root mean square 
(rms) dispersion is 0.13 M©, but SNe with intermediate decline rates [Ami5(5) ~ 1.05 to 1.5 magnitudes] 
show variations of almost a factor of 2 for the same value of Ami5(B). These SNe could cause scatter 
about the mean luminosity-decline-rate relation il6) . 

Once the contributions of ^^Ni, ^^Fe, and ^^Ni are added together to evaluate the total NSE mass, the 
dispersion decreases to 0.09 Mq. If SNe with different amounts of ^^Ni, and thus presumably different 
temperatures, but similar NSE content have similar AmisiB), it is likely that abundances (17) rather than 
temperature {18) primarily determine the opacity and light-curve shape. 

The outer velocity of the NSE region, determined from the width of the Fe lines, correlates with SN 
luminosity (I79l) . It coincides with the innermost silicon velocity, marking the transition from complete to 
incomplete burning: the difference between these two velocities, A{v) = 650 ± 900 km s"\ is consistent 
with zero. Remarkably, this results from two different methods applied to data obtained almost 1 year 
apart. 

Thus, while the mass of ^^Ni, and consequently the SN luminosity, can differ significantly, other 
characteristics of SNela are remarkably homogeneous. In particular, the narrow dispersion of the outer 
Si velocity indicates a similar extent of thermonuclear burning in all SNela: SNe that produce less ^^Ni 
synthesize more IMEs. 

The simplest interpretation of these seemingly antithetic results is that thermonuclear burning con- 
sumes similar masses in all SNela. We explore whether this scenario is consistent with the Mch model 
(liOl) . Applying the density-velocity structure of model W7 {Jj, we transform i'(Si) to mass (Fig. [T]). 
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Table 1 : Results and DDT parameters of parametrized three-dimensional delayed detonation models. 



model 




M(NSE) (Mo) 


M(IME) (Mo) 


toDT (s) 


Pddt (lO^gcm^^) 


D800 


1.004 


0.638 


0.547 


0.675 


2.40 


D20 


1.237 


0.833 


0.435 


0.724 


1.92 


D5 


1.524 


1.141 


0.220 


0.731 


1.33 



We find that the outer shell of siHcon encloses a mass of at least ~ 1.05 ± 0.09 Mo (~ 1.01 ± 0.05 Mo, 
excluding HVG SNe), independently of Am 15(B). This is a lower limit to the burned mass. 

The light-curve width t is related to /Sm\${B) (l20l) and depends on the ejected mass Mq, the kinetic 
energy E^, and the opacity /c as t oc k^/'^E'^^'^M^I'^ (|271) . We test our assumption computing parametrized 
light-curve widths and comparing them to observed values. The resulting luminosity-decline-rate relation 
(Fig. 121) is very tight, and practically identical to the observed one. This result supprots our opacity 
parametrization, corroborating our hypothesis that the mass burned is similar in all SNela. Given its 
weak dependence on E\^, the light-curve width is not much affected if more of the outer part of the WD 
is burned to IME, as these do not contribute much to the opacity. 

If SNela burn a similar mass, the progenitor mass is also likely to be the same, namely Mch- Because 
the outcome of the burning depends essentially on fuel density, a variation in iron-group elements pro- 
duction in Mch WDs requires different WD expansion histories. This in turn depends on the details of 
the burning. Once the WD reaches Mch, a thermonuclear flame is ignited near the center. The flame must 
start as a subsonic deflagration (i ), mediated by microphysical transport and accelerated by turbulence. 
As it propagates outwards, it could undergo a deflagration-to-detonation transition (DDT) and continue 
as a shock-driven supersonic detonation wave, in a so-called delayed detonation (1221) . This constitutes 
the most extreme explosion scenario admissible, exploring the limits of Mch explosions. 

We modeled the explosion using a three-dimensional level-set approach {23 ). The ignition of the defla- 
gration flame was treated as a stochastic process generating a number of ignition spots placed randomly 
and isotropically within 180 km of the WD centre (|2?1 |. 

What would cause the DDT is unclear. We assumed that it occurs as turbulence penetrates the internal 
flame structure: the onset of the so-called distributed burning regime (l25l) . This happens at low fuel 
densities, after some WD pre-expansion in the deflagration phase. The detonation is triggered artificially 
where the chosen DDT criterion is first satisfied, typically near the outer edges of the deflagration struc- 
ture (Table [J shows DDT parameters). Three simulations, with 800, 20, and 5 ignition spots, termed 
D800, D20, and D5, respectively, were performed on a moving cellular Cartesian grid (26) comprising 
the full star. Model D800, with its dense distribution of ignition points, exhausts the carbon-oxygen fuel 
at the WD centre almost completely in the deflagration phase. The energy release quickly expands the 
star, and the subsequent detonation mainly transforms low density outer material to IME. In contrast, 
the few ignition spots of model D5 consume little material during the deflagration, leaving more fuel at 
high densities which is converted mostly to NSE isotopes in the vigorous detonation phase. Model D20 
provides an intermediate case. 

The ejecta compositions of the model explosions agree grossly with the results derived from the spec- 
tra. The NSE mass produced ranges from 0.638 to 1. 141 Mq. The weakest explosions result from optimal 
burning in the deflagration phase (l26l) . and subluminous SNela are not reached in our parametrization. 
However, the conditions for DDT need further investigation. Eventually these events may be explained 



4 



within a single framework. Some extremely luminous SNela may come from very rapidly rotating WDs 
whose mass exceeds (27), but these are rare. 

The distribution derived from observations of burning products inside SNela could result from the 
variation of a single initial parameter, the flame ignition configuration, in Mch delayed detonations. 
The luminosity-decline-rate relation can be reproduced using this distribution and a simple opacity 
parametrization. Our results support the Mch scenario for most SNela, adding confidence to their use as 
distance indicators. 
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Figure 1 : The Zorro Diagram: Distribution of the principal isotopic groups in SNela. The enclosed mass (linked to velocity 
via the W7 explosion model) of different burning products is shown versus decline-rate parameter Ami5(S) (a 
proxy for SN luminosity). Individual SNe are coloured according to their velocity evolution {Hj: HVG, blue; low 
velocity gradient (LVG), green; and faint, red. Open circles indicate the mass of stable ^''Fe+^^Ni for each SN; 
solid circles that of '*Ni, and open triangles the sum of these (total NSE mass). Crosses show the sum of NSE 
and IME mass, indicating the total mass burned. The IME mass is the dilference between crosses and triangles. 
''*Fe and ''^Ni are found in roughly constant amounts in the deepest parts of all SNe, irrespective of luminosity: 
M(stableNSE) = 0.20 ± 0.05 Mq (lower horizontal shaded area). The ^*Ni mass determines the SN luminosity. It 
correlates with Ami5(6): M(^*Ni) = 1.34 - 0.67Ami5(B)[Me], with rms dispersion 0.13 Mq (lower diagonal shaded 
area). The total NSE mass correlates with hmii(B) better than M(^''Ni): M(NSE) = 1.55 - 0.69Ami5(B)[MG], rms 
dispersion 0.09 Mq (upper diagonal shaded area). IME lie mostly outside the iron-group zone. The outer Si velocity 
is similar for all SNe except HVG SNe. The mass enclosed by IMEs represents the total burned mass. When all 
SNe are included the average value is ~1.05 + 0.09 M© (upper horizontal shaded area). Excluding HVG SNe the 
value is ~1.01 ± 0.05 Mq (horizontal line). Both values are independent of iS.misi.B)- For a version of this plot with 
SN names and a bar diagram, see figs. S2 and S3. 
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Figure 2: Observed and synthetic luminosity-decline-rate relations for the SNe in our sample. Colors indicate velocity evo- 
lution as in Fig.[T] The peak luminosity L was computed from M(^''Ni) as L = 2 x lO'*^ M('*Ni) (|2S). Observed 
bolometric light-curve (LC) widths were obtained from observed rise and decline times (29) a&T = (t-i/i + f+i/2)- 
When /_i/2 is missing, it was estimated from the relation between ?_i/2 and Amis(B) derived from the other SNe. 
The dashed line is a linear regression between luminosity and observed LC width for the seven SNe common to our 
sample and that of i29\ . Synthetic bolometric LC widths were obtained assuming that: (i) The LC width t depends 
on ejected mass M^j, kinetic energy Zsk, and opacity k us: t <x k '-E^^^'^mI^'* (|23. (ii) For all SNe, Mej = Ma. The 
W7 GJ density-velocity distribution was used, (iii) The explosion kinetic energy depends on the burning product: 
Et = [1.56M('''^Ni) + 1.74M(stableNSE) + 1.24M(IME) - 0.46] x lO''' erg (77). (iv) Opacity is mostly due to line 
absorption i30\ . Accordingly, NSE elements contribute much more than IMEs because their atomic level struc- 
ture is more complex. The opacity was therefore parametrised according to the abundances of different species: 
K oc M(NSE) + 0. 1 M(IME) il6) . To compare the parametrized LC widths t to observed values T, a scale factor 
X = T It was computed for each of the seven SNe common to our sample and that of (29). The average factor, 
X = 24.447, was used to scale all SNe. Dots show the individual SNe. The continuous line is the linear regression 
between luminosity and synthetic LC width for our 23 SNela sample. For a version of this plot with SN names, see 
fig. S4. 
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